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Abstract 

We study thermodynamics of the apparent horizon in F(R) gravity. In particular, we demon- 
strate that a F(R) gravity model with realizing a crossing of the phantom divide can satisfy the 
second law of thermodynamics in the effective phantom phase as well as non-phantom one. 
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I. INTRODUCTION 



Recently, there have been more and more evidences to support that the current expan- 



sion of the universe is accelerating [l, 2]. The scenarios to explain the current acce 



to explain tne currer 
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erated 



121. One 



expansion of the universe fall into two broad categories 
is to introduce "dark energy" in the framework of general relativity The other is to study a 
modified gravitational theory, such as F(R) gravity, in which the action is described by an 
arbitrary function F(R) of the scalar curvature R (for reviews, see 
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n 



On the other hand, various observational data [13J imply that the ratio of the effective 
pressure to the effective energy of the universe, i.e., the effective equation of state (EoS) 
w c q = Pes/Pes, may evolve from larger than —1 (non-phantom phase) to less than — 1 
(phantom phase 14|). Namely, it crosses —1 (the phantom divide) at the present time or 
in near future. Recently, an explicit model of F(R) gravity with realizing a crossing of the 



151 ] . We note that the phantom crossing in 



phantom divide has been constructed in Ref. 
the framework of general relativity has also been studied in the literature, e.g., "quintom" 
model jic] ]. 

It is believed that a modified gravitational theory must pass cosmological bounds and 
solar system tests because it corresponds to an alternative theory of gravitation to general 



relativity. However, at the initial studies of F(R) gravity, models proposed in Refs. 



HQ, 



20| with the powers of the scalar curvature are strongly constrained. In recent years 



various investigations for viable models of F(R) gravity 
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331 ] have been executed. A curvature singularity problem in F(R) gravity has 



also been discussed in Refs. 



34 



As another touchstone of modified gravity, it is interesting to examine whether the sec- 
ond law of thermodynamics can be satisfied in the models of F(R) gravity. The connection 
between gravitation and thermodynamics was examined by following black hole thermody- 
namics (black hole entropy 36] and temperature [s3]) and its application to the cosmological 
event horizon of de Sitter space 38j. It was shown that Einstein equation is derived from 
the proportionality of the entropy to the horizon area together with the fundamental ther- 
modynamic relation, such as the Clausius relation 39(. This idea has been employed in a 



cosmological context 



411. 41 



42. 



431 ] . It was demonstrated that if the entropy of the apparent 



horizon in the Friedmann- Robertson- Walker (FRW) spacetime is proportional to the appar- 



2 



ent horizon area, Friedmann equations follow from the first law of thermodynamics 43]. 
The equivalent considerations for the FRW universe with the viscous fluid have also been 
studied 



44]. 



In addition, it was proposed [45| that in F(R) gravity, a non-equilibrium thermodynamic 
treatment should be required in order to derive the corresponding gravitational field equation 
by usin g th e procedure in Ref. (39)]. It was reconfirmed in Ref . [46 1 in F(R) gravity as well 



as Ref. 



471 ] in scalar-tensor theories. The first 48j and second 49j laws of thermodynamics 



on the apparent horizon in generalized theories of gravitation have recently been analyzed 
by taking into account the non-equilibrium thermodynamic treatment. Reinterpretations of 
the non-equilibrium correction 45j through the introduction of a mass-like function 50j and 



other approaches 



51 



52j ] have also been explored. Incidentally, the horizon entropy in four- 



dimensional modified gravity 53j and a quantum 



ogarithmic correction to the expression 



of the horizon entropy in a cosmological context |54j, [55|, [56| have been examined. Moreover, 
the first law of the ordinary equilibrium thermodynamics in F(R) gravity, scalar-tensor 
theories, the Gauss-Bonnet gravity and more general Lovelock gravity have been discussed 



in Refs. 



57 



58 



59 



60| , while the corresponding studies on the second law in the accelerating 



universe, F( R) g ravity, the Gauss-Bonnet gravity and the Lovelock gravity have been done in 



Refs. 



61 



62 



631 ] . respectively. Studies o: 



as well as its properties of dark energy 



thermodynamics in braneworld scenario 



64 



65 



66] 



671 ] have also been performed. 



It was pointed out in Ref. 68j that thermodynamics in the phantom phase usually leads to 
a negative entropy. Moreover, it was noted 69J that in the framework of general relativity the 
horizon entropy decreases in phantom models. However, the conditions that the black hole 



entropy can be positive in the F(R) gravity models 2l|, |22j, |23j with the solar-system tests 
have been analyzed in Ref. [70I ] . These recent studies have motivated us to explore whether 
in the framework of F(R) gravity the second law of thermodynamics can be satisfied in the 
phantom phase. To illustrate the point, in the present paper we consider a F(R) gravity 
model with a crossing of the phantom divide 15j as it contains an effective phantom phase. 
We use units of &b = c = U = 1 and denote the gravitational constant 8nG by n 2 = 8n/Mp 2 
with the Planck mass of M P1 = G~ 1/2 = 1.2 x 10 19 GeV. 

The paper is organized as follows. In Sec. II, we explain the first and second laws of 
thermodynamics in F(R) gravity. In Sec. Ill, we demonstrate that the model of F(R) gravity 
with the phantom crossing [ijj] can satisfy the generalized second law of the thermodynamics. 
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Finally, conclusions are given in Sec. IV. 



II. THERMODYNAMICS IN F(R) GRAVITY 

In this section, we study the first and second laws of thermodynamics of the apparent 
horizon in F(R) gravity. We consider the four-dimensional flat spacetime. 



A. F(R) gravity 

The action of F(R) gravity with matter is written as 

'F(R) 



d x\J—g 



2k 2 



(2-1) 



where g is the determinant of the metric tensor g^ v and £ ma tter is the matter Lagrangian. 
From the action in Eq. (12.11) . the field equation of modified gravity is given by 

F'(R)R, U - l -g^F{R) + g^UF\R) - V,V U F\R) = k 2 T^^ , (2.2) 

where the prime denotes differentiation with respect to R, V M is the covariant derivative 
operator associated with g^, □ = g^V^Vu is the covariant d'Alembertian for a scalar field, 
R^u is the Ricci curvature tensor, and j , j™ attcr ) — diag (p, p, p, p) is the contribution to the 
energy-momentum tensor from all ordinary matters with p and p being the energy density 
and pressure of all ordinary matters, respectively. 

We assume the flat Friedmann- Robert son- Walker (FRW) space-time with the metric, 

ds 2 = -dt 2 + a 2 (t)'y ij dx i dx i , (2.3) 
7ydxW = dr 2 + r W , (2.4) 

where a(t) is the scale factor and dVt 2 is the metric of two-dimensional sphere with unit 
radius. In the FRW background, from (//, v) = (0,0) and the trace part of (//, v) = 
(i, j = 1, • • • , 3) components in Eq. (12.21) . we obtain the gravitational field equations: 

H 2 = K ] z (P + Pc) , (2.5) 

K 2 

H = ^2F 7 {R) { ' P + P + Pc+Pc) ' (2 ' 6) 
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where p c and p c can be regarded as the energy density and pressure generated due to the 
difference of F(R) gravity from general relativity, given by 



Pc = — 

K z 



Pc = — 

K z 



- (-F(R) + RF'(R)) - 3HRF"(R) 
X - (F(R) - RF'(R)) + (2HR + F"(R) + R 2 F"'(R) 



(2.7) 
(2.8) 



respectively, with the scalar curvature of R — 6 + 2H 2 ^j. Here, H = a/a is the Hubble 
parameter and the dot denotes the time derivative of d/dt. We define the effective energy 
density and pressure of the universe as p c s = p t /F'(R) and p c s = p t /F'(R) with p t = 
p + p c and p t = p + Pc, respectively. Hence, from Eqs. ( I2.5jl and (12.61) we see that even 
in F(R) gravity, the gravitational field equations are expressed as H 2 = n 2 p c g/3 and H = 
—k 2 (p c fj +p c ff) /2, which are the same as those in general relativity. 

The continuity equation in terms of the effective energy density and pressure of the 
universe is given by 

p cS + 3H (p c{{ + p cS ) = 0. (2.9) 

Similarly, the (semi-)continuity equation of ordinary matters has the form 

p + 3H{p + p) =q. (2.10) 

One can take q = because the gravity is determined only by ordinary matters. Assuming 
that the energy fluid, generated from the modification of gravity, behaves as a perfect fluid, 
we have similar semi-continuity equations as 

p c + 3H (p c + p c ) = q c , (2.11) 
p t + 3H(p t +p t ) =q t , (2.12) 

where q c and q t (— q + q c ) are quantities of expressing energy exchange. Using Eqs. (12.51) . 
(ESD and (T2321 . we obtain 

3 2 dF'(R) 

<h--2 H —aT- (2-13) 

Clearly, from Eq. (I2.13p . we find that q t = for general relativity with F(R) = R, whereas 
qt does not generally vanish in F(R) gravity since there could exist some energy exchange 
with the horizon. 
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B. First law of thermodynamics 



We now illustrate the first law of thermodynamics in F(R) gravity. By using the spherical 
symmetry, the metric (I2.3P can be written as 



ds 2 = h a 3<ix a dx p + f W 



(2.14) 



where f = a(t)r, x° = t and x l = r, and h a @ is the two-dimensional metric h a p 
diag (—1, a 2 ). The dynamical apparent horizon is determined by the relation h a ^d a rdpf = 0. 
The radius of the apparent horizon for the FRW spacetime is given by 



13. 



i 

H 



(2.15) 



The associated temperature T of the apparent horizon is determined through the surface 
gravity of 

1 



AC 



d„ V^hh a(3 d, 



(2.16) 



2V-h 

where h is the determinant of the metric h a p. We note that the recent type la Supernovae 
data suggests that in the accelerating universe the enveloping surface should be the apparent 



horizon rather than the event one from the thermodynamic point of view 
In the FRW spacetime, one has 43] 



61] 



T 



AC. 



1 



2tt 2ixf A \ 2Hr A 
In general relativity, the entropy is expressed as 

A 



(2.17) 



S, 



GR 



AG 



(2.18) 



where A = 4wr\ is the horizon area. It follows from Eqs. (I2.17P and (12 . 1 8[) that 

TdSG * = h{ 1 ~ wrj dfA 

W dH 2 , W H 2 , 



(2.19) 



where V = 47rr^/3 is the volume of the apparent horizon. From Eq. (12. 5p . we see that H 2 = 
H 2 {p,p c ,F'{R)) and hence dH 2 /dt = {dH 2 /dp) p + {dH 2 /dp c ) p c + {dH 2 /dF'{R)) F'{R). 
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Substituting this equation into Eq. (12.191) and multiplying the resultant equation by the 
factor (k 2 /3) (dH 2 /dp)~\ we get 



TdS GR = -Vpdt - V 
dH 



dH' 



-1 #2 



-V 



dp 

- 1 ' dH 2 



dp J \dF'{R) 
The left-hand side of Eq. (I2.20p can be rewritten as 



F'(R)dt. 



dH' 
dp c 



p c dt 



k 2 f dH 2 
dp 



TdS< 



GR 



Td 



K 2 f dH 2 
dp 



GR 



dH' 
dp 



(2.20) 



(2.21) 



which leads to the Clausius relation 



TdS = 5Q , 



where the entropy S and the energy flux SQ are defined by 

k 2 fdH 2 ^ _1 



s 



3 V d P 
5Q = -Vpdt - V 



+ S GR d 



Sgr 

dH 2 
dp 



1 H 2 , fdH 2 
-dt-V 



H 



dH' 
dp 



2\ -1 



V 



dH 

dp 



dp 

2\ -1 



dH' 

dpc 

dH 2 
dF'(R) 



(2.22) 
(2.23) 

p c dt 

F'(R)dt. (2.24) 



From Eqs. (12.51) . (I2.18P and (12.231) . the entropy in F(R) gravity is expressed as 

AF'(R) 



(3s 



71 



72|: 



S 



4G 



(2.25) 



We note that there is arbitrariness in the definitions of S and 5Q. The reason for using 
those in Eqs. (12.231) and (I2.24p is that we can obtain the simple form of Eq. (12.251) as the 
expression of the entropy in F(R) gravity. Using p cS = 3H 2 /k 2 and Eq. (12. 5p . we find 



PeS — 9 



PcS + PeS 



dH' 

dp 

_j^ + A 

F'(R) k 2 
Pt + Pt 
F'(R) ' 



dH 2 \ . ( dH 2 

Pc + ( 1^777^ ) F (R) 



dpc 



dH' 
dpc 



\dF'(R) 
dH 2 



dF'(R) 



F'(R) 



(2.26) 
(2.27) 



where we have used (dH 2 /dp) = n 2 / (SF'(R)). By combining Eqs. (EM . fl2~TUD and QZZ 
we obtain 

.2 



dH 2 \ f dH 2 



dpc 



dF'(R) 



K 



3F'(R) 



[3H (p t + p t - p - p) + q] 



(2.28) 



Consequently, from Eqs. (I2l)]) . (l2~T0D . (l2~24l) and (T2~28D . we get 

5Q = 3H (p t + pt) Vdt - ^ (pt + Pt) ^ + T5 GR d (2.29) 
= -dE t + W t dV + ftWt + TSqrcI (F'(R)) , (2.30) 

where £? t = p t V is the total intrinsic energy and W t = — (1/2) T^ a ^h a /3 = (p t — pt) /2 
is the total work density {73]]. Here, = diag (pt,Pt,Pt,Pt) is the contribution to the 
energy-momentum tensor from all ordinary matters and energy fluid. This may be regarded 
as the work generated through the evolution of the apparent horizon [eJ. It follows from 
Eqs. (12.171) . ( 12.181) and (12.301) that the first law of thermodynamics in modified gravity can 
be constructed as 

5Q = -dE t + W t dV - Td c S , (2.31) 

where 

d c S = ~q t Vdt - S GR d (F'(R)) (2.32) 
8 ^ AH2 + H -d(F> m . (2.33) 



k2 (2H 2 + H^)H 2 

From Eq. (12.331) . it is clear that d c S does not vanish if F(R) is not equal to R. Hence, 
the emergence of the entropy production term is a special feature of thermodynamics in 
F(R) gravity. We note that such feature also appears in scalar-tensor theories because in 
those theories F'(R) is not constant but some dynamical quantities in terms of scalar fields. 
We remark that there exists a possible singularity in Eq. (12.331) if 2H 2 + H — 0, which 
corresponds to T = since T = (2H 2 + if) / {AttH) due to Eqs. (l2~T5l) and (|2~T7I) . It is 
clear that the necessary condition for a positive temperature is 2H 2 + H > 0. In this case, 
no singularity appears in Eq. ( 12.331 ). 

Using the Clausius relation in Eq. ( 12.221) . the first law of thermodynamics in F(R) gravity 
in Eq. ( 12.311) can be rewritten as 49] 



TdS + TdcS = -dE t + W t dV , (2.34) 

which characterizes the non-equilibrium thermodynamics of the apparent horizon in F(R) 
gravity. 
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C. Second law of thermodynamics 



Next, we investigate the second law of thermodynamics in F(R) gravity. From Eq. (12.341) . 
the first law of thermodynamics in terms of the horizon entropy Sh is expressed as 

TdS h = -dE t + W t dV - Td c S . (2.35) 

The Gibbs equation in terms of all matter and energy fluid is given by 

T t dS t = d (p t V) + Pt dV = Vd Pt + ( Pt + p t ) dV , (2.36) 

where T t and St denote the temperature and entropy of total energy inside the horizon, 
respectively. We assume that 

T t = bT, (2.37) 

where b is a constant with < b < 1. If there is no energy exchange between the outside 
and inside of the apparent horizon, i.e., q t = 0, thermal equilibrium realizes and therefore 
b= 1. 



The second law of thermodynamics in F(R) gravity can be described by 49] 



S h +^p- + S t >0 (2.38) 



[l-b)p t V + I l- h - ) ( Pt +Pt)V >0 (2.39) 



or 



2, 

by using Eqs. (l2T35l) - (l2~37l) . With Eqs. ([23]) and Q2~M . the relation in Eq. fl2~39l) is reduced 
to 

%h J *°> (2 ' 40) 

where 

J = (1 - b) H Z RF"(R) + 2(1-6) H 2 HF'(R) + (2 - b) H 2 F'(R) . (2.41) 

Thus, the condition to satisfy the second law of thermodynamics in F(R) gravity is equiva- 
lent to J > 0. 

We note that for F(R) = R, in which q t — 0, b — 1 and d c S = with the realization 
of thermal equilibrium, we find J = H 2 from Eq. (12.411) . In general relativity with a pure 
de Sitter expansion (H = 0), in which F(R) is given by F(R) = R — 2 A with A being the 
cosmological constant, J = 0. It is clear that J does not vanish for b — 1 if H ^ 0. 



III. SECOND LAW OF THERMODYNAMICS IN A F(R) GRAVITY MODEL RE- 
ALIZING A CROSSING OF THE PHANTOM DIVIDE 



In this section, we examine whether a F(R) gravity model [15( with the phantom crossing 



can satisfy the second law of thermodynamics discussed in Sec. II. In the model 
Hubble rate H(t) is given by 

27+1 ' 



15|, the 



7+(7+i) (y 



1 



27+1 



(3.1) 



where 7 is a positive constant and t s is the time when the Big Rip singularity appears. Here, 
we only consider the period < t < t s . When t — > 0, i.e., t <C t s , H(t) behaves as 

IO7 



H(t) 



t 



(3.2) 



In the FRW background Q, the effective EoS w cS is given by (8| w cS = -1 - 2H / (SH 2 ). 
In the limit of t — > 0, w c f[ = — 1 + 1/ (I57) > —1, corresponding to the non-phantom phase. 
The form of F(R) is given by 

^ v /60 7 (20 7 -l) J R- 1/2 "" 



F(R) 



1 



27+1 



R 



^607 (2O7- l)R~ l l 

57 ' Pj pj [6O7 (20 7 - l)] ft/2 R~^ /2 



><E 



2O7 



(3.3) 



where 



1 ± 



IOO7 (7 



(3.4) 



Here, to is the present time and p± are arbitrary constants. We remark that the stability 
for the obtained solutions in Eq. (13.31) u nder a quantum correction coming from conformal 



anomaly has been examined in Ref. 



151 ] . It has been shown that the quantum correction 



could be small when the phantom crossing occurs, although it becomes important near the 
Big Rip singularity. 



Using Eq. (JH2D, we find R ~ 6O7 (2O7 - 1) /t 2 and ^607 (2O7 - l)R~ l ' 2 /t s ~ t/t s < I. 
As the denominator inside the first large braces { } on the right-hand side of Eq. (I3.3P is 



10 



approximately unity, Eq. (13.31) can be simplified to 



F(R) 



-V6O7 (2O7 



->-> 



^-57/2+1 



5 7 -A 



2O7 



[6O7 (20 7 - l)] ft/2 iT^' /2 



(3.5) 



In this case, we obtain 
J = 5O7 2 



1^/607(207 - 1) 7 ^iT 5 ^ 2 ^{[10 7 (1-6) (-57 - ft + 2) - (2 - b)] 



x I ^q^' 1 ) ^ [607 (20 7 - l)] ft/a (5 7 + & - 2) iT^ 2 } . (3.6) 

Note that R > for # ~ 6O7 (2O7 - 1) /t 2 and 7 > 1/20. Using Eqs. fICTD . fl3"^|l and M) 
and taking into account the fact that p± are arbitrary constants, the necessary condition to 
have J > is 



2 - b\ > 



(3.7) 



P± {57 (1 - [(IO7 - 3) ± VI + IOO7 (7 + I) 

where we have assumed 7 > 1/20. We remark that for simplicity, we have chosen positive 
coefficients of i?~^ +//2 and R~^~l 2 in the large braces { } of Eq. (13.61) . By taking the values 
of p± so that the relation ( 13.71) can be met, the second law of thermodynamics, i.e., J > 0, 
can be satisfied. 

It follows from w e g = —1 — 2H / (3H 2 ) that w e s = — 1 when H = 0. Solving w e s = — 1 
with respect to t by using Eq. (13.11) . we find that w c s crosses the phantom divide at the time 
t = t c , given by 



-27+ a/47 



On the other hand, when t — » i s , we have 

H(t) 



1 



7 



1/(27+1) 



10 



t s -t 



(3.8) 



(3.9) 



In this case, the scale factor is given by a(t) ~ a (t s — t) 10 with a constant of a. When t — > t s , 
a — > 00 and therefore the Big Rip singularity appears. In this limit, w c ff = —16/15 < —1, 
corresponding to the phantom phase. The form of F(R) is given by 



F(R) 



( f 1 [t s -3VT^R^ 1/2 ] X ^ 



to 



1 _ 3vm R -i/2 



27+1 



Rj^Pj [t s - sVuor- 1/2 



ft 



X ti- 





ts + (P 3 - 15) R~ 1/2 



t s - 3VUOR-V 2 



(3.10) 
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which is reduced to 



F(R)~FR 7/2 (3.11) 



for t R ^> 1, where 



F =- 
7 



3^140(27 + 1) V*o 



t X7l5 



In this case, we obtain 



441000 (72 - 716) F 
(*, - «)" 

For 7 > 1/20, because i? > and < b < 1, the necessary condition to have J > in 
Eq. (I3TT3D is 

Efe^)>°> ( 3 - 14 ) 

which can be met by choosing the appropriate values of p-t. Since p± are arbitrary integration 
constants, for simplicity, we choose p + = 0. When 7 ~ 0(1) and p_ > 0, -F(-R) is always 
positive in both non-phantom and phantom phases. This is reasonable because for general 
relativity, F(R) = R > 0. Consequently, in this F(R) gravity model with the crossing of 
the phantom divide, the second law of thermodynamics can be satisfied. 

For a power-low type F(R) gravity described as F(R) = C\M 2 (R/M 2 ) n , where C\ and 
n are dimensionless constants and M denotes a mass scale, the scale factor a(t) is given by 
a(t) =a(t s - t )(™+ 1 )( 2 ™+ 1 )/(™+ 2 ) ^ ^j. T h e f orm Q f F ( R } in Eq . (pm]) corresponds to the 

case with n = —7/2. Accordingly, a(t) = a(t s — t)~ °, which implies a = 110a (t s — t)~ 12 > 
0. Thus, a late-time cosmic acceleration can be realized. This is the outcome of the F(R) 
gravity model in Eq. (13. lip . It should be noted that whether the gravity model of F(R) ~ 
FR 7 ' 2 in Eq. fl3~TT|) can pass Solar System tests still needs to be examined |75j |. 

From Eq. (13. 9p . we obtain R = 1260/ (t s — t) 2 , which leads to J = 
[(72 — 71b) F/4536] R 9 ^ 2 based on Eq. (13. 13D . In the thermal equilibrium limit, i.e., b ~ 1, 
we find J ~ (F/4536) R 9 ^ 2 . On the other hand, for general relativity with F(R) = R, it 
follows from Eq. (12.411) with 6=1 that J = H 2 = (1/15876) R 2 by assuming the same 
behavior of H in Eq. (13. 9p . The main difference between the expressions of J in the thermal 
equilibrium limit and that for general relativity is only the power of R. This comes from 

12 



the difference of the action between the present F(R) gravity F(R) ~ FR]I 2 in Eq. (13. lip 
and general relativity with F(R) = R. 

Finally, we remark that even in the effective phantom era of this F(R) gravity model, 
the second law of thermodynamics can be satisfied due to the non-equilibrium thermody- 
namic treatment. Hence, this model is more similar to a phantom model with ordinary 
thermodynamics suggested in Ref. [7^ . 

IV. CONCLUSION 

We have investigated the first and second laws of thermodynamics of the apparent hori- 



zon in F(R) gravity 49(. We have shown that in the F(R) gravity model with realizing 
the crossing of the phantom divide proposed in Ref. [ljj], the second law of thermodynamics 
can be satisfied in not only the non-phantom phase but also the effective phantom one. In 
addition to cosmological constraints and solar system tests on the models of F(R) gravity, 
such an examination whether the second law of thermodynamics can be met in those models 
is important. The demonstration in this work can be regarded as a meaningful step to con- 
struct a more realistic model of F(R) gravity, which could correctly describe the expansion 
history of the universe. 
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